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ABSTRACT 



Aims. Using XMM-Newton data, we have aimed to study the nuclear outflow of the nearby starburst galaxy NGC 253 in X-rays with 
respect to its morphology and to spectral variations along the outflow. 

Methods. We analysed XMM-Newton RGS spectra, RGS brightness profiles in cross-dispersion direction, narrow band RGS and EPIC 
images and EPIC PN brightness profiles of the nuclear region and of the outflow of NGC 253. 

Results. We detect a diversity of emission lines along the outflow of NGC 253. This includes the He-like ions of Si, Mg, Ne and O 
and their corresponding ions in the next higher ionisation state. Additionally transitions from FeXVII and FeXVlII are prominent. 
The derived temperatures from line ratios along the outflow range from 0.21 ± 0.01 to 0.79 ± 0.06 keV and the ratio of FeXVIl lines 
indicates a predominantly collisionally ionised plasma. Additionally we see indications of a recombining or underionized plasma in 
the FeXVII line ratio. Derived electron densities are 0.106 + 0.018 cm"^ for the nuclear region and 0.025 ±0.003 cm"^ for the outflow 
region closest to the centre. The RGS image in the O VIII line energy clearly shows the morphology of an outflow extending out to 
~750 pc along the south-east minor axis, while the north-west part of the outflow is not seen in O VIII due to the heavy absorption 
by the galactic disc. This is the first time that the hot wind fluid has been detected directly. The limb brightening seen in Chandra 
and XMM-Newton EPIC observations is only seen in the energy range containing the Fe XVII lines (550-750 eV). In all other energy 
ranges between 400 and 2000 eV no clear evidence of limb brightening could be detected. 

Key words. X-rays: galaxies - Galaxies: individual: NGC 253 - Galaxies: spiral - Galaxies: starburst - ISM: jets and outflows 



c3 . 1. Introduction 



Starburst galaxies are known to show very complex emission in 
X-rays. This emission originates on the one hand from sources 
that appear to be point-like sources, like X-ray binaries, su- 
pernovae and supernova remnants. On the other hand emis- 
sion comes from the diffuse hot component of the interstel- 
lar medium like diffuse emission in the disc and gaseous out- 
flows driven out of the disc by massive stellar winds and core 
collapse supernovae, also called superwinds. T he latter phe- 
nome non can be quite spectacular i n M82 (e.g . IStevens et alj 
12003b and NGC 253 (e.g. iStricklandetalJ|200a), for example, 
where the se superwinds emerge from a starburst nucleus. In 
NGC 253. lFabbiano & TrinchiCTild 19841) first detected the south- 
east part of this outflow in X-rays with E instein and called it 
the "minor-axis component". With j?(9&4r. |Pietsch et"an ( |2000|) 
also detected the part of the outflow pointing into the opposite 
direction. However the spatial resolution of ROSAT was not yet 
good enough to learn more about the morph ology of this outflow . 
Later on, observations with X MM-Newton dPietsch et alJl200ll) 
and especially with Chandra dStrickland et al. " l2000h showed 
that the outflow can be explained with a Umb brightened hollow 
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cone structure. Temperatures of the best fit thin thermal plasma 
models are in the range 0.15-0.94 k eV from XMM-Newton EPIC 
and 0.46-0.66 keV from Chandra. [Strickland eTaTI (|2000|) con- 
cluded that the detected emission originates from the shocked re- 
gion at the border of the outflow where the wind collides with in- 
terstellar medium. The wind itself though was thought to be too 
hot and too thin to be detected directly. This picture, however, 
disregards the possibility that the wind may be mass-loaded, en- 
training ambient interstellar medium (ISM) as well as infalling 
material. If turbulent mixing proceeds on a time scale which is 
larger than the flow time within a given region, such as the base 
of the outflow studied here, we expect some dumpiness in the 
outflow, imprinted on an overall less dense wind. As we will 
show later, this is confirmed by our analysis, which shows that 
e.g. O VIII is not limb brightened. 

High resolution spectra of NGC 253 and M82, taken with the 
XMM-Newton Reflection G r ating Spect rometer (RGS) were firs t 
published by iPietsch et all (l2001h and iRead & Steveiisl ( |2002|) . 
respectively. Both spectra show the LyQ. emission lines from Si, 
Mg, Ne, O, N and also their helium-like charge states. Also 
both galaxies show emission lines from Fe XVII and Fe XVIII, 
and M82 additionally shows lines from Fe XX, Fe XXIII and 
Fe XXIV. In M82 the line ratios for neon, iron and oxygen are 
quite dilFerent compared to NGC 253. In general the M82 spec- 
trum appears to be hotter with temperatures in the range of ~0.3- 
1 .5 keV, with its continuum more confined to higher energies. Its 
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X-ray flux, as well as its X-ray luminosity, in the RGS energy 
band (0.35-2.5 keV) is higher than that of NGC 253. 

However these spectra only give a combined spectrum of the 
nuclear source and the outflow. In this paper we present an anal- 
ysis where we decompose the total spectrum of NGC 253 into 
regions containing the nucleus and different parts of the outflow, 
while maintaining the high spectroscopic resolution. 



2. Observations and data reduction 

The nuc lear region of NG C 253 was observed with XMM- 
Newton jJansen et al.l 1200 ll) during two orbits in June 2000 
and June 2003 using all of the European Photon Imaging 
Camera (EPIC) instruments ( IStriider et alj 120011 : [Turner et al.l 
1200 lb and the two co-aligned RGS spectrometers, RGSl and 
RGS2 (Iden Herder et al.ll200lh . for a total of about 190 ks. The 
details of these observations are shown in Table[T] An additional 
archival observation (Obs. id. 0110900101) could not be used 
for this analysis, since the pointing of this observation was in 
the north-west halo of NGC 253 with the result that the outflow 
was not in the field of view (FOV) of the RGS. 

Before we start to describe the analysis procedures, we want 
to place emphasis on why it is possible at all to perform the 
following spectroscopic analysis. First, the nuclear outflow of 
NGC 253 is an extended X-ray object, which can be spatially 
resolved by XMM-Newton EPIC and RGS, embedded in an even 
larger region of X-ray emission from point-like sources and dif- 
fuse emission in the disc and halo of the galaxy. This means the 
RGS data for the central regions are in principle affected by the 
contamination from the surrounding emission, but, as shown in 
FigU both the nucleus and the outflow in NGC 253 are signif- 
icantly brighter and well above the galaxy emission, so we can 
expect only a minor contamination. Moreover as we show later, 
we can identify and "remove" effects due to the disc emission. 

Second, since the RGS is a slitless spectrometer, the spectra 
of all sources in the field of view are superimposed on each other 
on the detector. A spatial displacement of a source along the dis- 
persion direction corresponds to a wavelength shift in the spec- 
trum of 2.31 X 10"^^ A arcsec with respect to a not-displaced 
source. Since the outflow has an extent of up to 1 .4' in the dis- 
persion direction, the spectral resolution is limited to ~0.19 A at 
15 A. This is still considerably better than the energy resolution 
from CCD detectors. 

Due to the superposition of all sources, other bright point 
sources in the FOV could contaminate the spectrum of the out- 
flow. However the effective area decreases significantly for off- 
axis sources and even a contribution from t he bri g htest off-axis 
source at the botto m of Fig. |2](X21 from iPietsch et al.i 120001 : 
iTanaka etalj|2005[) can be neglected. The source at the south- 
west edge of the outflow in region SE 1 (cf. Fig.[T|and|2]| however 
does affect the outflow spectrum. T he spectrum of this source 
does not show line features (X33 in lPietsch et al.ll200lh . so its 
contribution to the RGS spectrum from this region is an in- 
creased continuum flux. This does not affect our conclusions. 

We analysed the data using the Science Analysis System 
(SAS), version 6.5.0, together with the most recent calibration 
files available at the time of the analysis. The metatask rgsproc 
1.19.6 was used to process the RGS data. We first filtered the 
data sets for time periods with low contamination by low energy 
protons. Then we extracted a light curve from the background 
region on CCD 9, the chip closest to the optical axis of the tele- 
scope, therefore the most affected by background flares, to de- 
termine the threshold count rate, which we then used to filter 
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Fig. 1. EPIC PN brightness profiles along the RGS dispersion di- 
rection of each extraction region. The strong peak at 0" is caused 
by the outflow emission. In the Centre region, shown in Fig.|2] 
this is superimposed by the nuclear source of the galaxy. Point 
sources in the extraction regions are seen as sharp spikes in the 
profile. A positive distance points parallel to the major axis to- 
wards the south-west. The distance from the galactic major axis 
is given by the value z in the captions of the individual extraction 
regions. 



the eventfiles. The thresholds were 0.25 ct s ' for observations 
0125960101 and 0125960201, and 0.20 ct s"' for observation 
0152020101, where the count rate was more stable during the 
non-flaring time intervals. 



2.1. RGS spectra 

Since the dispersion direction of the RGS was aligned approxi- 
mately with the major axis of the galaxy we were able to extract 
spectra for different adjacent regions along the minor axis, i.e. 
the cross-dispersion direction, of the galaxy (see Fig.|2l). The ex- 
tent of the extraction regions are 30" in cross-disper sion direc- 
tion. A ssuming a distance to NGC 253 of 2.58 Mpc (iPuche et al.l 
Il99ll) this corresponds to a width of the extraction regions of 
375 pc. The events in these regions were additionally filtered 
with a CCD pulseheight filter to select only the m - —\ spectral 
order The reference points for the origins of the energy scales of 
the spectra were set on the minor axis of the galaxy, where the 
outflow has its peak emission (see Fig.[TJ. 

NGC 253 is an extended source and covers most of the area 
of the RGS detectors. To prevent contamination in the back- 
ground spectra we used the task rgsbkgmodel 1.1.5, which 
computes background spectra from RGS background templates. 
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Table 1. XMM-Newton NGC 253 observation log. Observation identification, observing date, pointings and orientation of the satel- 
lite, total exposure time (Texp) and exposure time after screening for high background (Texp, dean) are given. 



Nr. 


Obs. id. 


Obs. dates 


Pointing 


direction 


P. A. 


T 

^ exp 


T 

-'- exp, clean 








RA/DEC (J2000) 


(deg) 


(ks) 


(ks) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


1 


0125960101 


2000-06-03 


00:47:36.74 


-25:17:49.2 


56.9 


60.8 


45.1 


2 


0125960201 


2000-06-04 


00:47:36.57 


-25:17:48.7 


57.0 


17.5 


7.0 


3 


0152020101 


2003-06-19 


00:47:36.89 


-25:17:57.3 


53.8 


113.0 


75.9 




Fig. 2. Logarithmically scaled EPIC PN-hMOS image of the central region of NGC 253 in the energy band 0.5 - 2.0 keV. The centre 
of the galaxy is marked with the black cross. The D25 ellipse is overlaid in black, brightness contours and the extraction regions for 
the RGS spectra are overlaid in white. The black X at the bottom of the image marks the source X21. 



To increase statistics we combined the spectra of the two 
RGS detectors and added up the spectra from the three obser- 
vations. As the position angles of XMM-Newton did not differ 
very much in the three observations, the regions from which the 
spectra were extracted are only slightly tilted relative to each 
other A difference in position angle between observations has 
the effect of degrading the spatial coincidence in the extraction 
regions. However in our case position angle differences are small 
enough that we can neglect this error. A correction of the spec- 
tra for effective area and the combination of the spectra from 
different observations and instruments was done with the task 



rgsfluxer. The task's description states that the fluxed spec- 
trum produced by rgsfluxer should not be used for any serious 
analysis of the data. If we use it in spite of this warning, we have 
to consider the following effect that will add to uncertainties: 
The task rgsfluxer neglects the redistribution of monochro- 
matic response into the dispersion channels, so the intrinsic line 
broadening of the detector is not removed from the spectrum. 
As we do not determine line positions or line widths and since 
we integrate over the whole line including its wings to derive 
fluxes for individual lines, this effect does not restrict our anal- 
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ysis. Additionally we are only interested in relative line fluxes 
and not in absolute values. 

To obtain acceptable statistics (> 3cr) for most of the 
emission lines while maintaining a high spectral resolution of 
~0.39 A in the spectra we combined six channels at a time. 

For the observations the dispersion direction of the spec- 
trometers was approximately aligned along the major axis of the 
galaxy. Therefore the spectra for different cross-dispersion areas 
correspond to regions with different distance z from the galactic 
plane (cf. Fig. |2]i. The spectra obtained from these regions are 
shown in Fig.|3] 

2.2. RGS images 

The RGS is a slitless spectrometer So the whole observed target 
within the FOV is imaged on the detector plane, with an offset in 
dispersion direction for every wavelength it emits in. This aspect 
can be used to extract narrow band images for various emission 
lines. In cross-dispersion direction the image in a selected emis- 
sion line is directly mapped onto the detector CCD, while in dis- 
persion direction the observed object is compressed into a nar- 
row region. The technique was fir st applied to RGS data of th e 
supernova remnant DEM L 7 1 by Ivan der Hevden etai] (|2003|) . 
With the help of K. J. van der Hey den (private communications) 
we developed our own code to produce these narrow band im- 
ages. The procedure is as follows: The low background eventlists 
are filtered for the wavelength range of the desired line and the 
"banana region" in wavelength-energy-space to exclude second 
order spectra and noise. By setting the wavelength range narrow 
enough we made sure that no neighbouring lines would be in the 
extraction region. The dataset thus obtained has to be converted 
into spatial coordinates and it has to be uncompressed along the 
dispersion axis using the following equation as described in the 
SAS task rgsangles: 



Afi. 



sin a F 

A0- 

siny6 L 



(1) 



with the change in the grating exit angle Afi due to the offset 
in the angular component A(f> of an off-axis source parallel to 
the dispersion direction, the angle of incidence a, the grating 
exit angle /?, the focal length F, and the distance between the 
Reflection Grating Assembly and the prime focus L. 

The images were corrected for exposure and binned to a 
pixel size of 0.4". In a next step we included a RA-DEC co- 
ordinate system. As reference coordinate, we chose the coordi- 
nate of the centre of NGC 253 (the position of the brightest IR 
source in the galaxy, Qf2ooo = 0^7™ 33^3, ^2000 = -25°17'18", 
iForbes et al . 2000). In cross-dispersion direction the position of 
our reference coordinate on the CCD could be taken directly 
from the source list file which was produced in the processing of 
the RGS data. In dispersion direction the position is given in the 
above procedure by the Doppler shift corrected line centre po- 
sition. The velocity that had to be accounted for is the ga laxy's 
systemic velocity of 243 km s ' dKoribalski et al.ll2004l) . This 
shifts the reference coordinate along the dispersion direction by 
~0.35" A"' times the centre wavelength in which the image is 
calculated, e.g. 6.6" for the O VIII image. The effects of the ve- 
locities of the earth with respect to the sun, and XMM-Newton's 
orbital velocity can be neglected since they are only of the order 
of 0.7" and 4"xl0""^, respectively in the OVIII image where 
the effect would be largest, and therefore much smaller than the 
width of the point spread function. The images for the lines were 
created separately for each of the three observations and were 



then combined into one image. In a final step the images were 
smoothed with a gaussian filter. 

In general the method is affected by two different effects: (i) 
A Doppler shift due to the radial velocity component of an emit- 
ting source changes its position in the image along the dispersion 
direction axis. A radial velocity of 1000 km s ' would corre- 
spond to 21.6"at a wavelength of 15.0 A. (ii) Assuming there is 
not only one but two lines in the wavelength extraction interval, 
we would have two images of the object in the resulting image, 
superimposed with an offset along the dispersion direction of 
7.2' A"'. For example the emission in the two FeXVII lines at 
16.780 A and 17.055 A would be superimposed with an offset of 
~2.0'. The Fe XVII at 17 A image is the only case where we ac- 
tually have to consider that we have created an image using two 
lines, i.e. the FeXVII lines at 17.055 A and at 17.100 A. In the 
spectra the lines are not separated and they appear about equally 
strong, which is mos t likely an effect d ue to the low statistics. 
According to theory (iMewe et al.ll 1985b the line strength of the 
latter line should be ~52-85 % of the former one, depending on 
the temperature of the plasma (1.08 keV and 0.11 keV, respec- 
tively). The separation of the lines is 0.046 A which corresponds 
to a shift in dispersion direction of 19.9". By smoothing the im- 
age with a larger gaussian with a FWHM of 20" we can account 
for the error we make by using both lines. 

Images in the NeX, FeXVII and OVIII lines are shown 
in Fig. |5] The number of photons that were extracted from all 
three observations are 806, 1077, 816 and 1231 for the images 
in the Ne X (1 1.98-12.35 A), FeXVII at 15 A (14.86-15.13 A), 
FeXVII at 17 A (16.90-17.21 A) and OVIII (18.80-19.17 A) 
lines respectively. The FWHM of the gaussian filter was 12" for 
Ne X, Fe XVII at 15 A, O VIII, and 20" for Fe XVII at 17 A. 



2.3. RGS cross-dispersion profiies 

To get additional spatial information of the line distribution, we 
produced emission line profiles in the cross-dispersion direction. 
Therefore we extracted events from the RGS eventfiles by apply- 
ing the same filters in wavelength and wavelength-energy-space 
as for generating the RGS images, but then we binned the counts 
into 30" bins to match the extraction regions that were applied 
for the spectra. Background counts in the respective wavelength 
ranges were taken from the spectra that were obtained with the 
RGS background model task and subtracted from the emission 
line profiles. Four of these profiles are shown in Fig.|4] 



2.4. EPIC-PN images 

To verify the results from the RGS images we also extracted 
EPIC PN narrow band images in approximately the same en- 
ergy ranges. Therefore we filtered the PN eventfiles in the en- 
ergy bands around Ne X (992-1052 eV), Fe XVII at 15 A (795- 
844 eV), OVIII (625-690 eV) and FeXVII at 17 A (694- 
734 eV). The spectral resolution of the EPIC PN detector is 
~70 eV, so it is possible that photons with higher or lower ener- 
gies contribute to the energy band of interest. Additionally there 
is contamination from higher energies due to the redistribution 
in the detector. Photons can lose up to 60% of their energy in the 
CCD before they are detected. This means that bright features in 
some energy range can show up to some degree in lower energy 
bands. The filtered eventfiles of the different observations and 
instruments were merged using the SAS task merge. We created 
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images of the eventfiles and smoothed them with a gaussian of 
6". The resulting images are shown in Fig.|6] 

2.5. EPIC-PN brightness profiles 

To detect the Umb br ighte ning of the outflow , as found by 
[Strickland eTaTI (|2000|) and iPietsch et all (1200 ll) . we also ex- 
tracted brightness profiles perpendicular to the outflowing di- 
rection from the merged PN eventfiles of all observations. 
Furthermore to check for an energy dependence in the limb 
brightening we subdivided them into energy bins with a width of 
~ 150 eV starting from 400 eV up to 2000 eV. The emission lines 
of O VIII and the Fe XVII lines are included in the energy ranges 
550-700 eV and 700-850 eV, respectively. Ne X is mostly in the 
1000- 11 50 eV bin. The extraction regions match the regions we 
used for the RGS spectra. Furthermore we split the region SE 1 
into the two regions 'SE 1 (1)' and 'SE 1 (2)' with a width of 15" 
in cross-dispersion direction of the RGS. The brightness profiles, 
sorted for energy band, are shown in Fig.|7] Since we are not in- 
terested in the emission from the nucleus nor from the bright 
source X33 south-west of the nucleus, the profiles are limited to 
a maximum of 250 counts. This still shows the main features in 
the Centre and SE 1 (1) region, while cutting off the peaks in 
some cases. 



3. Results 

3.1. RGS spectra 

The RGS spectra show emission in many different lines (Fig.|3]l. 
Especially in the region including the major axis of NGC 253 
(Centre) a large variety of lines from different elements can 
be identified. The spectra extend from the Si XIV line at the 
highest energy down to the C VI line at the low energy end 
(A = 6... 34 A). All the He-like ions in this range (Si XIII, 
Mg XI, Ne IX, O VII and N VI) and their corresponding ions in 
the next higher ionisation state (Si XIV, MgXII, NeX, OVIII 
and NVII) can be resolved. The iron 3d-2p transitions around 
15 A as well as the 3s-2p lines around 17 A are detected, i.e. 
their peak heights are more than twice the error in the wave- 
length bin. Also the lines at ~16.0 and ~16.1 A can be clearly 
detected, but it is not clear whether the line at ~16.0 A is from 
Fe XVIII or from O VIII. In the range 10 A to 1 1 .5 A there is an 
indication of iron lines from Fe XXIII and Fe XXIV. 

A spectral feature that is prominent in the spectra is a broad 
base at the FeXVII lines at 15 A and 17 A and at the OVIII 
and O VII line positions with a width of up to ~0.5 A. If this line 
broadening is caused by the velocity dispersion of the outflowing 
gas, it would imply deprojected velocities of up to ~40000 km/s, 
corresponding to temperatures of ~ 10" K. Such high veloci- 
ties are in disagreement with outflow velocities in other galax- 
ies. Extreme cas es show values of up to ~3600 km/s (e.g. in 
NGC 3079, Veill eux et al.|[T994.) . A more plausible explanation 
for the broad base is contributions from the disc emission of 
NGC 253: In Fig.[T]disc emission extends about 200" in disper- 
sion direction (major axis of the galaxy) in both directions from 
the centre. This corresponds to a ~0.5 A shift in wavelength and 
therefore can explain the observed effect. 

The region north-west of the centre (NW) is strongly affected 
by absorption from the galactic disc that lies between the out- 
flow and the observer .Pietsch et al. (2000) derive an additional 
absorbing column A^h of 1 - 2 x 10^^ cm"^ for this position 
north-west of the centre. Therefore most of the lines are only 



weak or not detectable at all. If we assume that the north-west 
and south-east outflow have a similar intrinsic spectrum, then 
absorption can fully account for the difference between regions 
NW and SE 1, at the same projected distance from the centre. 
The strongest lines in the NW region are from NeX, NelX, 
FeXVII at -15 A and -17 A, OVIII and from the forbidden 
line of the O VII triplet, but most of the emission originates from 
the disc and is smoothed out into the broad base, as mentioned 
above. 

The regions south-east (SE 1 to SE 3) of the centre do not 
suffer from this absorption by the disc. One can easily follow 
how the lines increase or decrease in strength as one is going 
away from the galactic disc south-east along the minor axis. The 
lines at short wavelengths from silicon and magnesium are the 
first to disappear with distance from the centre. The Ne X line 
is seen to decrease considerably in strength south-east of the 
Centre region, whereas the line from the lower ionised NelX 
is not affected as much and is even stronger than the Ne X line 
in regions SE 1 and SE 2. All the lines from iron decrease in 
strength except for the FeXVII lines at 17 A which grow by a 
factor of -1.5. Also the OVIII line increases in strength. The 
O VII triplet has about the same strength in region SE 1 as in 
the Centre region. The lines from N VII and C VI increases in 
strength compared to the central region. 

Further away from the centre, in region SE 2, the FeXVII 
lines at ~15 A are still detectable and the OVIII line is the 
strongest line in the spectrum. Also the OVII triplet is still 
strong. The detection of all the other lines is below 2cr, even 
though the lines at wavelengths longer than 12 A can still be 
identified. 

In region SE 3 only very weak lines from O VIII and from 
the Ne IX triplet remain. 

Unfortunately the statistics in the spectra are not good 
enough to allow a quantitative spectral analysis with XSPEC. 
However, several conclusions can be drawn. 

Temperature estimates can be inferred from line ratios of dif- 
ferent elements or of the same element in different ionisation 
states using model calculations. Assuming collisional ionisation 
equilibrium (CIE). lMewe et al] (1 1985b calculated line strengths 
for different elements and transitions depending on the temper- 
atures of the plasma. By measuring the fluxes of two transitions 
in a spectrum and taking the ratio of these, the obtained value 
can be compared with the tables in iMewe et"al] (1 19851) and the 
temperature of the plasma can be derived. 

We used the line strength ratio between the Ly,, state of a 
given element and its helium-like charge state that matched the 
morphology of the Ly„ state in the RGS images best. The line 
strengths were derived by integrating the flux of the line over the 
wavelength (cf. Table |2]i. We found that the derived temperature 
value and its error show only a weak dependence on how much 
of the wings of the line we include in the flux integration. The 
resulting temperatures for Si, Mg, Ne and O and their variation 
along the outflow direction are shown in Table [3] Using only 
the peak height of the line however gives temperatures that are 
lower, except for oxygen in the regions Centre and SE 1, com- 
pared to the values shown in Table[3]by up to 50%. 

Using the line flux from an emission line (cf. Table|2]i and an 
estimate of the size of the emitting region we can derive electron 
densitie s for the nuc l ear re gion of NGC 253 and the south-east 
outflow. iMewe et"ar l (fT985h give the line power P' normalized to 
the electron density for different temperatures and X-ray emis- 
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Energy [keV] 

2.066 1.536 1.222 1.015 0.868 0.758 0.675 0.605 0.549 0.505 0.464 0.451 0.402 0.577 0.554 
0.6 




6 7 8 9 10 1112 15 14 15 16 17 18 19 20 21 22 25 24 25 26 27 28 29 50 51 52 55 54 55 

Wavelength [A] 

Fig. 3. Combined RGS spectra of NGC 253 extracted from different regions along the outflow. The label gives the south-east position 
z of the extraction region along the minor axis relative to the centre of the galaxy in arcmin. 
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Ta ble 2. Flux values f or d ifferent emission line s in the extraction regions. The references for the expected wavelengths (/lexpected) 
are lMewe et al] ( Il985h and lPhilUps eTal] ( Il999l) . 
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3.0 ± 2.0 
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13 ±3.9 




OVIII 


19.0 


6.7 ± 2.5 


16 ±3.8 


18 ±4.0 


6.9 ± 2.5 


OVII w 


21.6 




2.3 ± 2.2 


3.0 ± 2.4 




OVII x+y 


21.8 




3.4 ± 2.6 






OVIIz 


22.1 


5.1 ±3.3 


9.5 ± 4.4 


13 ±5.2 




NVII 


24.8 




2.2 ± 1.5 


6.3 ± 2.8 


1.5 ± 1.3 


NVI w 


28.8 


1.4 ± 1.3 




4.2 ± 2.4 




NVI x+y 


29.1 






2.4 ± 1.8 




NVIz 


29.5 






3.2 ± 2.0 


1.8 ± 1.5 


CVI 


33.7 


3.0 + 2.5 


4.5 ± 3.0 


6.6 ± 3.3 


5.7 ± 3.1 



Table 3. Temperatures of the plasma for different regions of the 
outflow of NGC 253 derived from line ratios of different ele- 
ments. For region SE 3 no temperatures could be derived as the 
lines are too weak. 



Region 



Temperature in keV 





Si Mg 


Ne 


O 


NW 


0.61±0.08 


0.51 ±0.08 


0.21±0.01 


Centre 


0.79±0.06 0.66±0.04 


0.43±0.02 


0.22±0.01 


SE 1 


0.46±0.04 


0.38±0.03 


0.21±0.01 


SE2 




0.25±0.02 


0.31±0.04 



sion lines. The electron density can then be derived using the 
formula 



iFEyAncP- 



VP' 



(2) 



where iig is the electron density, F the flux (in counts s"' cm"-) 
in an emission line with Ey from an emitting region with volume 
V and distance d. The emitting region is assumed to be uniform 
in «„. 



For the central region we selected the Ne X line and assumed 
a uniform, spherical emitting region with a radius of ~160 pc. 
The outflow in region SE 1 is best represented in the OVIII 
line and we assumed a uniform, cylindrical volume with a ra- 
dius of ~200 pc and a height of 375 pc, the latter being con- 
fined by the extraction region. We selected these lines because 
they are strong in the spectra and because we can get a good 
estimate for the emitting volume from the RGS images. The re- 
sulting electron densities are «e,nucieus - 0.106 + 0.018 cm"^ and 
"e.outflow = 0.025 + 0.003 cm"-' for the nucleus and the outflow 
region, respectively. 

Apart from the derivation of temperatures and electron den- 
sities, some selected emission lines can be used as diagnostic 
lines. 

The FeXVII lines at 15 A and 17 A can be used to de- 
rive the ionising mechanism in the plasma. In the Centre region 
of NGC 253 as well as in the regions NW and SE 2 the line 
strengths indicate a predominantly collisional ionised plasma. 
Region SE 1, however, shows an inverted line ratio. Here the 
lines at 17 A are stronger than the lines at 15 A, which points at 
a photoionised plasma. 

In general the helium-like line triplet of O VII can provide 
the electron de nsity, the electron te mperature as well as the ioni- 
sation process ( iPorquet et al.l200TI) . However, we refrained from 
using the OVII triplet for the following reasons: In the com- 
bined and fluxed spectra the significance is below 3 cr for most of 
the spectral bins. Also the individual lines are not clearly distin- 



8 



M. Bauer et al.: High-resolution X-ray spectroscopy and imaging of the nuclear outflow of the starburst galaxy NGC 253 




-50 50 1 00 150 -50 50 100 150 

Distance [orcsec] Distonce [arcsec] 

Fig. 4. Ne X (upper left), Fe XVII at 15 A (upper right), O VIII (lower left) and Fe XVII at 17 A (lower right) profiles of background- 
subtracted counts against cross -dispersion distance for the combined NGC 253 RGSl and RGS2 data. The bin at distance zero 
corresponds to the extraction region 'Centre'. Negative distances are towards the north-west, and positive values correspond to 
areas to the south-east. 



guishable from each other. The latter is probably enhanced due 
to the use of the task rgsfluxer as described above. We also re- 
frained from doing a simultaneous fit of the single uncombined 
spectra with XSPEC, as the statistics in one single spectrum are 
barely above 2 cr for the strongest bin. 

3.2. RGS cross-dispersion profiles 

In the cross-dispersion profiles (Fig.|4|i line emission is strongest 
in the Centre bin at a distance of 0" for NeX and FeXVII at 
15 A. In the O VIII and FeXVII at 17 A profile this is not the 
case. Here the SB 1 region at a distance of +30" is the bright- 
est. However the north-west half of the central bin is already 
affected by absorption from the disc. When one bins the bright- 
ness profile into smaller spatial regions, the strength drops dra- 
matically in the central bin and towards the north-west. Towards 
the south-east, on the other hand, it only drops slowly. Further 
to the south-east with distances larger than 45" the Ne X and the 
FeXVII profiles show a large drop in brightness, whereas the 
O VIII profile declines with a less steep gradient. Therefore the 
Ne X and Fe XVII emissions are more concentrated within the 
disc, whereas the O VIII emission extends farther away from the 
disc. The general behaviour that the emission from higher en- 
ergy lin es is more concentrated and not as extended is also seen 
in M82 dRead & Stevensll2002l) . 

3.3. RGS images 

The RGS images allow us to derive additional information about 
the spatial composition of the outflow. But before we describe 
the images in more detail we need to point out that these images 



are dominated by low statistics. We therefore restrict our results 
and conclusions on the regions with good statistics. 

The strong OVIII emission clearly traces the outflow. It 
reaches out to ~750 pc projected distance along the south-east 
minor axis and has an extent perpendicular to that of ~400 pc. 
Because of the high absorption, emission from the farther side 
of the disc is not detected. With the RGS spatial resolution, the 
outflow in O VIII does not show limb brightening. The emission 
is strongest close to its central axis and its intensity decreases 
towards the border. This suggests that the emitting ions are not 
concentrated on the borders of the outflow where the outflow gas 
interacts with the surrounding material, but that this emission 
comes directly from the outflowing gas. The extent to the south- 
east is less than the one in the EPIC images where the bright 
outflow emission can be followed to ~1.3 kpc (Fig.|6]l, due both 
to the sensitivity and to the smaller energy band (AA =1.87 A 
vs. 0.37 A) that was used to extract the images for the RGS. 
The emission from Ne X is clearly concentrated in the disc and 
nucleus and it gives no contribution to the outflow. The images 
in the iron lines both show the strongest emission south-east of 
the nucleus and the emission is slightly extended along the out- 
flow direction. However they do not trace the outflow morphol- 
ogy as seen in the OVIII line. This is not surprising, because 
the excitation cross sections for the iron and O VIII lines have a 
different temperature dependence. Due to the work that is per- 
formed when the outflow expands against the pressure of the 
ambient medium, as well as to the divergence of the flow in the 
cone perpendicular to the disc, the temperature, or more pre- 
cisely the kinetic energy of the ele ctrons, should decrease (e.g. 
iBreitschwerdt & Schmutzlerll 19991) with height z above the disc, 
explaining the relative increase in emission of O VIII to iron. 
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O VIII 
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Fig. 5. RGS images of NGC 253 in the Ne X, Fe XVII and O VIII lines. The white lines mark the major (north-east to south-west) 
and minor axis (north-west to south-east) of NGC 253. The image in the Fe XVII line at 17 A was smoothed with FWHM of 20", 
all others with 12". White c ontours indicate 2cr, 3cr, 4cr, . . . above the background. The white cross south-west of the nucleus marks 
the bright source X33 from lPietsch et alj (l2000h . The RGS dispersion direction is such that wavelength increases from south-west 
to north-east as indicated by the arrow. 



3.4. EPIC-PN images 

We find that the EPIC FN narrowband images (Fig.|6]l are all af- 
fected by the redistribution effect of the detector, i.e. the image 
is contaminated by events with higher energies from outside the 
energy filter boundaries. Especially very bright sources like the 
central source of the galaxy contribute strongly to this effect and 
all images will therefore show these sources. The bright nuclear 
source is clearly visible in all of the EPIC FN images. Also the 
bright source ~0.5' south-west of the nucleus is clearly visible in 
all EFIC FN images. For the RGS images the latter is only true 
for the FeXVII at 15 A image. This clearly shows the advan- 
tage, namely the far better energy resolution, of the RGS images 
compared to the EFIC FN narrowband images. 

With respect to the limb brightening of the outflow emis- 
sion, the image in the energy range 694-734 eV (including the 
FeXVII at 17 A line) indicates a morphology that could result 
from a limb brightened outflow at a distance of ~0.7' away from 
the galactic centre. Also in the energy range 625-690 eV (includ- 
ing the O VIII line) there is an indication of this morphology at 
a distance of about 1.1'. The corresponding RGS image unfor- 



tunately has too low statistics to confirm this. The other images 
give no indication of limb brightening. 



3.5. EPIC-PN brightness profiles 

In the EFIC-FN brightness profiles (Fig. |7]i the peak in region 
SE 1 (1) that is caused by the outflow emission coincides always 
with the centre. A double peak with a central depression is visi- 
ble in SE 1 (2) and SE 2 at 700-1000 eV (covering also the iron 
lines). This may point at a limb brightened outflow in this region. 
An indication of this structure in region SE 2 at 550-850 eV 
could be a redistribution effect of the detector (see above). At 
other energies no indication of limb brightening can be seen. 
The profiles show either a flat plateau or a hump peaking around 
the centre of the outflow. This confirms the findings of the RGS 
images. Starting from 700 eV region NW comes out and it is 
stronger than the two SE 1 regions together in energies above 
1150 eV. The strong peak in the region SE 1 (1) (blue, dotted 
line) at ~+20" distance is the bright point source south-west of 
the nucleus and not part of the outflow region. 
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Fig. 6. EPIC PN images of NGC 253 in the energy bands around the Ne X (922-1052 eV), Fe XVII (795-844 eV and 694-734 eV) 
and O VIII (625-690 eV) lines. The images show the region of NGC 253 that is covered by the RGS images (Fig. |5]l. The thin 
black lines mark the major (north-east to south-west) and minor axis (north-west to south-east) of NGC 253. The area within the 
thick black lines gives the extraction region in the corresponding RGS line image. Black contours indicate 2cr, 3cr, Act, . . . above the 
background. The images were smoothed using a gaussian with a FWHM of 6". 



To rule out that our findings are affected by XMM-Newton's 
spatial resolution we compared brightness profiles, extracted 
with the same regions and energy bands from Chandra obser- 
vation 3931, to the EPIC PN profiles. There the double peak is 
clearly detected in region SE 1 (2) between 700 and 1000 eV, 
possibly also up to 1 150 eV. It is not visible in regions SE 1 (1) 
or SE 2. The profiles in region NW show the same behaviour 
as in the EPIC-PN data. So taking the differences in the instru- 
ments into account XMM-A^ewfon and Chandra give a consistent 
picture. 



4. Discussion 

4.1. Line ratios and temperatures 

The extracted RGS spectra of the outflow along the minor axis of 
the galaxy show emission lines from many ions in different ion- 
isation states: the Ly^ lines from Si, Mg, Ne, O, N and also their 
helium-like charge states. Additionally we see emission lines 
from Fe XVII and Fe XVIII. 

With increasing distance from the nuclear region, the rela- 
tive flux in the OVIII line intensifies compared to the flux in 
FeXVII, i.e. the line flux at longer wavelengths increases. The 
same effect can also be seen in the line flux ratio of O VII to 
O VIII. For the oxygen line ratio this implies that the temper- 
ature decreases; the gas is cooling as it flows away from the 



nucleus. This change in temperature would also affect the ex- 
citation of Fe XVII. Since both elements, oxygen and iron, are 
from the same kind of sources, i.e. type II SNe, it seems unlikely 
that the distribution is different, unless there is a very different 
dumpiness in O and Fe initially. Therefore a change in the oxy- 
gen to iron line flux is more likely due to the change in temper- 
ature than due to a different radial abundance profile. 

The line strengths of the Fe XVII lines at 15 A compared to 
the one at 17 A indicates a predominantly coUisional ionised 
plasma. Region SE 1, however, shows an inverted line ratio. 
There, the flux in the FeXVII lines at 17 A is enhanced com- 
pared to the lines at 15 A. This points at a photoionised plasma, 
but no strong photoionising sources were detected in the vicinity. 
There are two additional alternatives to create this line ratio: 

A) the plasma is highly underionised compared to the ionising 
electrons. The time since the heating of the plasma was too short 
to reach an equilibrium state. In these so called und erionised 
plasmas inner-shell ionisation is highly operational jKosenkd 
20061) and leads to an enhanced 17.10 A flux from FeXVII 
( Doron & BehaHl2002l) . Examples that s how underionised p las- 
mas are su pernova remnants like N 1 32 D jBehar et al.ll200ll) and 
DemL71 Jvan der Hevden et alj|2003h . 

B) the plasma is overionised. It can be produced in fast adiabatic 
cooling of hot {T ~ 10** K) and al most completely ionised gas 
expan ding out of a superbubble dBreitschwerdt & Schmutzlen 
Il999h . The ionisation cross sections as well as the recombina- 
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Fig. 7. EPIC PN brightness profiles of the outflow region for the RGS extraction regions sorted by energy band. Note that region 
SE 1 was split into the two regions 'SE 1 (1)' and 'SE 1 (2)' with a width of 15" in cross-dispersion direction of the RGS. The 
abscissa gives the distance from the minor axis of the galaxy in arcsec. Positive values are to the south-west, negative values to 
the north-east. The presented histograms consist of raw detector counts, i.e. the background and the detector characteristics, like 
quantum eflSciency or effective area, were not taken into account. 



tion cross sections are different for all ions, so when a shock 
propagates through the outflow region, the ionisation of dif- 
ferent elements and different ionisation stages occurs at differ- 
ent time scales for each. As the expansion of the wind occurs 
on much shorter time scales than radiative recombination, the 
highly ionised atomic states remain frozen-in. These states then 
recombine at a later time further away from the nucleus. This is 
called delayed recombination. In both regions. Centre and SE 1, 
we see emission from Fe XVIII, which is produced in the hot 



starburst region of the galaxy. As the wind breaks out of the nu- 
clear superbubble and expands into the lower density medium 
along the minor axis of the galaxy it cools mostly adiabatically. 
However the Fe XVIII in the wind recombines on a larger time 
scale than the wind fluid, so we have Fe XVIII in abundance 
although the electron temperature would preferentially admit 
FeXVII. When the Fe XVIII ions then recombine to FeXVII, 
the states of FeXVII which decay via the 17 A lines are prefer- 
ably populated compared to the ones which decay via the 15 A 
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lines. This enhances the flux at 17 A and could lead to the mis- 
leading line ratio in this region. 

Temperatures derived from line ratios of individual elements 
are in the range of 0.21 + 0.01 t o 0.79 + 0.06 keV. This range 
is larger than the one found by IStrickland et" ar (20QCk_ usmg 
Chand ra (0.46!°- - 0.66!jJi° keV). However StricklandetaD 
(I2OOOI) used only a single temperature MEKAL hot plasma 
model for each extraction region. Given the complexity of the 
X-ray producing mechanisms this approach is definitely too sim- 
ple. Our der i ved te mperature range is closer to that given by 
iPietsch etaLl(l2001l) fmmXMM-Newton EPIC spectra, who used 
three MEKAL models (0.15, 0.53 and 0.94 keV). The approach 
in this paper to derive temperatures tries to account for the high 
complexity of the outflow, without the claim to be complete, by 
using actually one MEKAL model for each element, and by tak- 
ing advantage of the superior energy resolution of the RGS. This 
approach gives already almost the w hole measured energy r ange 
just for the Centre region, whereas [Strickland et al.l (l2000l) ob- 
tain this range of temperatures by using several different extrac- 
tion regions along the outflow. Looking at the temperatures of 
just one extraction region we see that different elements give dif- 
ferent temperatures. This could be due to sampling of different 
regions along the line of sight, but also due to delayed recom- 
bination. The latter is a true non-equilibrium ionisation (NEl) 
situation and the result would be measuring different tempera- 
tures for different ions. All of the above temperatures were de- 
rived by the use of models that assume CIE for the emitting gas. 
This is not necessarily the case for an outflow. To avoid this as- 
sum ption we would have to fit the spec tra to NEl models (e.g. 
from lBreitschwerdt & Schmutzleilll999l) . These models entail a 
higher level of complexity as they depend crucially on a detailed 
hydrodynamical model of the outflow. Therefore they are sen- 
sitive to the time evolution of the wind, and as such unique. In 
other words, NEl models would enable us to fully exploit the 
information buried in the X-ray spectra, since one would calcu- 
late a series of models with different boundary conditions, and 
then derive suitably binned synthetic spectra for fitting the obser- 
vational results. Thus a satisfactory fit would determine within 
certain limits the physical parameters of the outflow. 

The Fe lines could in principle be used as well to derive 
temperatures. However the emission from FeXVll at ~15A 
(2p-3d) and at ~17A (2p-3s) are affected by delayed recom- 
bination of Fe XVIII, which contradicts the CIE assumption. 
Therefore the results could be misleading and are not shown 
in Table [3] However from the fact that the FeXVII at 15 A 
lines are present and strong in the spectra, the temperature has 
to be above 0.22 keV in the regions where the Fe L shell 
originates from, otherwise Fe XVII line formation is inefficient 
jRugge & McKenzielll985]) . Furthermore calculations show that 
above 0.60 keV Fe XVIII is strong while Fe XVII is weak. This 
is not the case in any of the regions where we obtained spectra 
from, therefore we expect that kT=0.2 . . . 0.6 keV 

The RGS images give another indication of the cause of the 
large temperature spread. The different distribution of the events 
for each emission line strongly points to the fact that plasmas 
in different regions of space are responsible for the line emis- 
sion. It is not surprising that these plasmas then do not have the 
same temperature and that the temperatures for each element are 
different. The time scale for attaining pressure and hence tem- 
perature equilibrium in the starburst region is of the order of the 
sound crossing time. If we compare this to the average time in- 
terval between supernova explosions, we can assess the smooth- 
ness of the temperature distribution. A rate of 0.05 yr"' has been 



quoted for NGC 253 dColina & Perez-01eal[T992l) . A supernova 
blast wave expanding within a hot tenuous gas has a low Mach 
number, so that the sound crossing time, Tsc, in the starburst re- 
gion is roughly the time scale for the shock wave to reheat the 
hot gas. Adopting a value of 20" or d = 250 pc and a temper- 
ature of 0.6 keV or 6.6 x 10^ K for the central region yields a 
value of Tsc ~ d/cs ^ 8.2 x 10^ yr, which is much larger than 
the interval of 20 years between successive explosions, even if 
we decrease the supernova rate by a factor of 10 - 100 for the 
starburst region. Therefore temperature inhomogeneities in the 
starburst region and hence variations in the ionisation stages are 
to be expected. As the flow moves away from the disc sources, 
the temperature structure should however become smoother with 
time when the flow time becomes larger than the sound crossing 
time Tsc- 

IStrickland et"an ( |2000|) argue that the spectral variation in 
Chandra spectra 'along and between the northern and southern 
outflow cones is due to variations in the absorption column and 
not due to significant temperature variations along the outflow'. 
A higher absorption column is certainly the case for the NW re- 
gion, resulting in too high temperature values. Although from the 
analysis of the XMM-Newton RGS spectra we cannot draw con- 
clusions regarding the absorbing column, we find that there is a 
significant temperature variation along the outflow. The temper- 
atures for Mg and Ne are decreasing with distance to the south- 
east from the centre. The temperature of oxygen seems to be 
constant, except for the value in region SE 2. 

4.2. The morphology of the outflow 

The only RGS image that clearly shows the outflow geometry 
is the one in the OVIII line. The shape of the south-east out- 
flow matches roughly the truncated cone with a n opening angle 
of ~2 6° that is seen in Chandra and H,j data dStrickland et al.l 
l2000h . though the angle in the OVIII image appears to be 
slightly smaller. Also this emission is not limb brightened. The 
image suggests that the outflow cone is filled with clumpy dis- 
tributed O VIII. If this is correct, then we see the emission from 
the hot wind fluid itself for the first time. The dumpiness may 
be caused by mass-loading, i.e. the turbulent process of mix- 
ing in ambient ISM and infalling material. This detection of 
the wind fluid however is in contrast to the statement from 
IStrickland & StevensI (l2000l) that the wind is too thin to emit ef- 
ficiently enough to be detected. 

In XMM-Newton EPIC obse rvations bvlPietsch et al.l (1200 ih 
and in Chandra observations bv IStrickland et al.l ( l2000l) the out- 
flow shows a limb brightened morphology. In our analysis we 
find limb brightening in the FeXVII at 17 A RGS images. The 
comparison between the different, though with lower spectral 
resolution, PN images and brightness profiles also shows an en- 
ergy dependence of limb brightening and a clear detection of the 
emission from the outflowing gas at the energies below 700 eV. 

IStrickland & StevensI (I2OOOI) found from simulations that 
low opening angles for outflows are only possible when a thick 
galactic disc is present. This produces opening angles from a few 
degrees when the wind starts to blow out of its superbubble and 
then increases up to values of ~ 60°, depending on boundary 
conditions. Thin-disc models typically show opening angles of 
~ 90°. However these simulations assumed an ISM distribution 
in rotating hyd rostatic e quilibrium with the gravitational field. 
Simulations bv Ide AviTl ez & Breitschwerdt (2004, 2005) show 
that the ISM in the disc has a highly complex structure. This af- 
fects the break-out dynamics of a superbubble via density and 
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pressure gradients and also via "holes" from previous outflows 
that have not closed yet. Conclusions about the thickness of the 
disc and its collimating effect are hence difficult. 

4.3. Outlook 

In a forthcoming paper we will analyse the XMM-Newton EPIC 
data of the outflow and the halo emissio n with the NEI model 
from iBreitschwerdt & Schmutzleij (Il999h . which will then also 
consider the ionisation history of the outflowing gas from the 
nucleus into the halo. This will give further information on the 
composition, temperature and dynamical structure of the out- 
flow. 
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